Quantum information approach to the Ising model: Entanglement in chains of qubits 
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Simple physical interactions between spin- 1/2 particles may result in quantum states that exhibit 
exotic correlations that are difficult to find if one simply explores state spaces of multi-partite 
systems. In particular, we present a detailed investigation of the well known Ising model of a chain 
(ring) of spin-1/2 particles (qubits) in a transverse magnetic field. We present explicit expressions 
for eigenstates of the model Hamiltonian for arbitrary number of spin-1/2 particles in the chain in 
the standard (computer) basis and we investigate quantum entanglement between individual qubits. 
We analyse bi-partite as well as multi-partite entanglement in the ground state of the model. In 
particular, we show that bi-partite entanglement between pairs of qubits of the Ising chain (measured 
in term of a concurrence) as a function of the parameter A has a maximum around the point A = 1 
and it monotonically decreases for large values of A. We prove that in the limit A — > oo this state is 
locally unitary equivalent to an JV-partite Greenberger-Horn-Zeilinger state. We also analyse a very 
specific eigenstate of the Ising Hamiltonian with a zero eigenenergy (we denote this eigenstate as 
the X-state). This X-state exhibits the "eXtreme" entanglement in a sense that an arbitrary subset 
A of k < n qubits in the Ising chain composed of N = 2n + 1 qubits is maximally entangled with 
the remaining qubits (set B) in the chain. In addition we prove that by performing local operation 
just on the subset B one can transform the X-state into a direct product of k singlets shared by 
the parties A and B. This property of the X-state can be utilised for new secure multi-partite 
communication protocols. 

PACS numbers: PACS numbers: 03.67.-a, 03.65.Ud,05.50.+q 



I. INTRODUCTION 



Those multi-partite quantum systems which are funda- 
mental objects of statistical and solid state physics, have 
been found interesting also from a perspective of quan- 
tum information processing. These systems often ex- 
hibit multi-partite entanglement that can be used cither 
for quantum information processing or quantum com- 
munication. Amongst such systems a distinguished role 
is played by exactly solvable models, such as the Ising 
model describing a chain of interacting spin-1/2 parti- 
cles in an external magnetic field. Eigenstates of the 
corresponding model Hamiltonian can be studied from a 
perspective of quantum information theory with a good 
physical motivation: Any quantum computer is a phys- 
ical device composed of elementary units, qubits, de- 
scribed by a certain Hamiltonian. Consequently, perfect 
knowledge of the Hamiltonians and their eigenvectors are 
vital. An important condition the physical system has 
to fulfil is the possibility of preparation of an a priori 
known initial state. The easiest way to realize this task 
is to simply let the system evolve into its ground state. 
Thus the knowledge of the entanglement properties of 
the ground state or more practically thermal states are 
necessary. This has been followed by many authors. In 
particular, various versions of the Heisenberg model (XX, 
XY, XYZ) have been investigated. Many of these studies 
concern numerical and analytical investigations primar- 
ily focused on the behaviour of bipartite entanglement 
of small number of qubits in ground and thermal states, 

e.g. Refs. fliiaiiaasiiiamiii The no - 

tion of "thermal entanglement" , i.e. the entanglement of 



thermal states is introduced, and its properties including 
threshold temperatures and magnetic field dependence 
are studied. 

Symmetry properties of multi-partite systems have 
been used to calculate entanglement among their con- 
stituents. In Ref. (13] thermal equilibrium states of 
isotropic two-spin systems are analysed exploiting SU(2) 
invariance. The results are related to isotropic Heisen- 
berg models. In Ref. 0], analytical expressions for cer- 
tain entanglement measures are derived using general 
symmetries of the quantum spin system. Then they are 
used for the XXZ model in order to calculate concurrence 
and the critical temperature for disentanglement for fi- 
nite systems with up to six qubits. It should be noted, 
that they use the 3-tangle to analyse some multi-partite 
entanglement aspects of the system, and discuss entan- 
glement sharing in detail. 

In Ref. the authors pointed out, that in a finite 
chain of qubits, the time evolution generated by the 
Ising Hamiltonian produces "entanglement oscillations" , 
which lead to the presence of GHZ and W type entan- 
gled states. A generalisation to 2D and 3D models is 
also outlined. Discussions of multi-partite entanglement 
also appear in Refs. In Ref 01 quantum tele- 

portation is utilised as a tool to reveal the importance 
of multi-qubit entanglement in a 3 qubit Heisenberg- XX 
chain. 

A central question in the problem of entanglement of 
more than two systems is that of bounds on entangle- 
ment. Three or more quantum systems cannot be arbi- 
trarily entangled in the similar way as they cannot be 
arbitrarily classically correlated ^21 • The state with an 
a priori specified entanglement properties may not exits 
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at all 0, an d therefore the search for a state with 
given, in a sense optimal entan glem ent properties, is in 
general a hard problem. In Refs. |2fil2ll | the authors have 
solved such particular problems by minimising the energy 
of a Hamiltonian. That is, the sought state with a given 
pairwise entanglement is the ground state of a Hamil- 
tonian with a very clear physical interpretation. Since 
such states may be useful for quantum information pro- 
cessing it is desirable to known the concrete form of the 
states that are either optimal, or obey certain bounds. 
It should be pointed out that the problem of finding en- 
tangled webs with give n properties has been extensively 
addressed in Refs. |18L l22j| but without any reference to 
systems described with Hamiltonians. Thus it is inter- 
esting to see how this issue can be approached in other 
exactly solvable models. 

Another interesting issue concerning such models of 
many-body systems is the collective behaviour of these 
systems under certain conditions known as the critical 
phenomena. Let us point out that these phenomena has 
already been studied extensively. On the other hand it 
has been pointed out only recently |23ll24ll25ll26ll27l | that 
entanglement is the quantity that may play a crucial role 
in the description and understanding of critical phenom- 
ena. The central concept of the theory of critical phe- 
nomena is the universality - the critical exponents char- 
acterising divergences near critical points are the same 
for all systems belonging to the same universality class. 
For a special class of one-dimensional magnetic systems 
it has been shown in Ref. |2^| that the bipartite entangle- 
ment shows scaling behaviour near the transition point. 
One should also expect, that precursors of the critical 
behaviour may appear even in non-critical systems. 

In this paper we present a detailed investigation of 
the well known Ising model of a chain (ring) of spin-1/2 
particles (qubits) in a transverse magnetic field (Section 
II). We present explicit expressions for eigenstates of the 
model Hamiltonian for arbitrary number of spin-1/2 par- 
ticles in the chain in the standard (computer) basis and 
we investigate quantum entanglement between individual 
qubits (Sections III and IV). We analyse bi-partite as well 
as multi-partite entanglement in the ground state of the 
model. In particular, we show that bi-partite entangle- 
ment between pairs of qubits of the Ising chain (measured 
in term of a concurrence defined in Section II) as a func- 
tion of the parameter A has a maximum around the point 
A = 1. In addition, it monotonically decreases for large 
values of A. We prove that in the limit A — > oo this state 
is locally unitary equivalent to an iV-partite Greenberger- 
Horn-Zeilinger state (Section IV). We also analyse a very 
specific eigenstate of the Ising Hamiltonian with a zero 
eigenenergy (we denote this eigenstate as the V-state). 
This X-state exhibits eXtreme entanglement in a sense 
that a arbitrary subset A of k < n qubits in the Ising 
chain composed of N = 2n + 1 qubits is maximally en- 
tangled with the remaining qubits (set B) in the chain. 
In addition we prove that by performing local operation 
just on the subset B one can transform the X state into 



a direct product of k singlets shared by the parties A and 
B. This property of the X state can be utilised for new 
secure multi-partite communication protocols. Techni- 
cal details of some of our calculations are presented in 
appendices. 

II. SETTING-UP THE SCENE 

A. The Ising model 

We consider a model of a linear chain of spin one-half 
particles forming a circle, placed in a magnetic field where 
only the z-component of the field is non-zero. Since we 
are interested in the spin degrees of freedom only, the 
Hamiltonian of the system is given by 

JV JV 
n— 1 n— 1 

where cr",o; = x,y,z are well known Pauli operators. 
The first term in the Hamiltonian is the interaction term 
with coupling constant Ci and the second term corre- 
sponds to a free Hamiltonian. The lower index n labels 
the position of a spin in the chain and N is the overall 
number of particles. The cyclic boundary conditions 

a 'N+i= a i> a = x,y,z (2.2) 

ensure that the chain forms a circle. The form of the 
interaction is chosen such that each particle interacts only 
with its two nearest neighbours. 

The Hamiltonian in Eq. (|2.1|) can be rewritten into a 
form which is more convenient for numerical calculations: 

{N N ~\ 

n=l n=l ) 

where E = B and A = Cj /B is a dimensionlcss parame- 
ter. Now, instead of two parameters in energy units (B 
and Cj) we have one dimensionless parameter A and one 
parameter in energy units, E, that can be neglected in 
our further calculations except for the investigation of 
entanglement in Gibbs states (see Section ITlI G|) . 

B. Measures of entanglement 

In this paper we will use three different measures - 
the concurrence, the tangle and a measure of an intrinsic 
three-partite entanglement . 

The concurrence |28j is a measure of the bipartite en- 
tanglement between two qubits. Let pab be the joint 
density matrix of the system consisting of qubits A and 
B. The matrix pab Pab has four non- negative eigen- 
values {Ai, A2, A3, A4} that are written in a descend- 
ing order (i.e. {Ai > A2 > A3 > A4}. The matrix 
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Pab is a spin-flipped version of density matrix pab-, i-e. 
Pab = o y ® (J yP*AB a y ® a v The concurrence is given by 
the relation 



Cab = max • 



;0 



(2.4) 



Let us point out that the state pab is arbitrary, that is, 
the concurrence is a valid measure of entanglement for 
two-qubit mixed states as well. The minimal value of the 
concurrence is zero (in this case two qubit states are pre- 
pared in a separable state), while for any Cab > two 
qubits are entangled. The maximal value of Cab = 1 
corresponds to maximally entangled states that arc lo- 
cally unitary equivalent to Bell states. It has been shown 
earlier that the concurrence is directly related to the en- 
tanglement of formation |2^| . 

On the other hand the tangle has been originally de- 
fined for pure states only. (There is an extension to mixed 
states but the extension is not computationally feasible 
except for the case of two qubits when the tangle is equal 
to the square of the concurrence). Consider a multi- 
partite system where one of the subsystems, denoted A, 
is a qubit. The tangle C a a between a subsystem A and 
the rest of the system, denoted as A reads 

C AA =4Det M = 2(l-Tr^), 



where pa is the reduced density operator describing a 
state of the subsystem A. 

Finally, the intrinsic three-qubit entanglement is de- 
fined for pure states only. Consider a system composed 
of three qubits A, B and C and let the system of three 
qubits be in a pure state. With the help of the tangle and 
concurrence introduced above we can define pure three- 
partite entanglement 



Cabc = C A a ~ C A b ~ C A c • ( 2 - 6 ) 

In Ref. 01 it has been proven that the definition (|2.6(l 
of intrinsic three-qubit entanglement is independent of 
permutations of particles and shares all properties that 
a proper measure of entanglement has to fulfil. 

Finally, let us mention that measures of entanglement 
are not unique and different measures might result in dif- 
ferent ordering of states. For the case of bipartite systems 
AB prepared in a pure state the measure of entanglement 
is in fact any suitable function of the eigenvalues of the 
reduced density matrix of either of the two subsystems 
Pa or pb Ref. [29|. For example the well known Von 
Neumann entropy 



S(pa) = -Tr p A log pa 



(2.7) 



defines a bipartite measure of entanglement usable for 
arbitrarily dimensional systems. 



III. THREE SPIN-1/2 PARTICLES 

In order to understand entanglement properties of the 
Ising chain under consideration it makes sense to start 



with a relatively simple example of three spin-1/2 parti- 
cles (qubits). Even this simple system exhibits interest- 
ing properties and their understanding will guide us in 
general case of an arbitrary number of qubits. 

The Ising Hamiltonian with three spin-1/2 particles 
can be directly diagonalized and energy levels easily cal- 
culated. In what follows, we will call spin-1/2 particles 
as qubits since the Hilbert space Ti. of a spin one-half 
particle is two-dimensional. Let us note that the sim- 
plest example is the case of two qubits, that is N = 2. 
However, this trivial example has already been investi- 
gated in Refs. and [j| . The case of three qubits is also 
interesting on account of the fact that, besides intrinsic 
bipartite entanglement, three qubits can also share three- 
partite entanglement. In the case of three qubits being in 
a pure state this intrinsic three-qubit entanglement can 
be easily calculated with the help of Eq. I|2.6|l . Finally, 
even such a simple example nicely illuminates main re- 
sults concerning multi-partite entanglement where most 
results can be generalised to the case with an arbitrary 
number N of qubits (spin-I/2 particles) in the chain. 



The Hamiltonian Eq. (|2.3[1 of the Ising model with only 
three qubits in the chain reads 



(2 5) H 3 = _ ^°f ® a 2 - Aof ® of - \af (g) of + a\ + + of . 

Note that the Hamiltonian H3 is permutationally invari- 
ant, unlike the Hamiltonians Hm for N > 3 qubits. All 
Hamiltonians if/v are obviously translationally invariant. 
What is not so obvious is the fact that the Hamiltonian 
is invariant under the inversion of the order of particles. 
The particles in the chain are labelled with n = 1, . . . N. 
Now, if we relabel them as n - -> A — n + 1 the Hamil- 
tonian remains unchanged and thus is invariant under 
the inversion of the order. For the case of N = 3 the 
two transformations, translation and inversion of the or- 
der, together with an arbitrary combination of the two 
yield in fact all possible permutations of the particles in 
the chain. Thus it follows that the Hamiltonian H3 is 
permutationally invariant. 

The knowledge of the symmetry of the Hamiltonian is 
utmost important, as it plays a crucial role in the pro- 
cess of finding its eigenvalues and eigenvectors. It is well 
known that for every symmetry S there exists a unitary 
or an anti-unitary operator T$ such that the correspond- 
ing Hamiltonian commutes with Ts 



[Ts,H}=0 



(3.2) 



As a result of this commutation relation the two oper- 
ators H and Ts have common set of eigenvectors. It 
means, that there is one set of vectors (basis of the cor- 
responding Hilbert space) which are eigenvectors of H as 
well as the operator T<j. Moreover, any non-degenerate 
eigenstate of H has to be invariant under the action of the 
operator T5. On the other hand, any eigenstate which is 
not invariant under the action of the operator Ts is de- 
generate. In what follows the knowledge of symmetries 
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FIG. 1: The spectrum of the Hamiltonian H3. We present 
the dependence of eigenenergies as a function of the coupling 
parameter A. The energy levels E\ = E2 and £3 = E4, are 
degenerate. The other four energy levels are non-degenerate. 
The ground state corresponding to the state with the lowest 
energy in our notation is represented by the seventh level E7 
for all values of the parameter A. 

of Hn will help us to find some particularly interesting 
states of the spin-chain under consideration. 

A. Spectrum of the Hamiltonian 

The spectrum of the Hamiltonian H 3 can be easily 
calculated directly by diagonalising the Hamiltonian H3. 
The Hilbert space H.2®'H.2®'^i-2 of three qubits is eight di- 
mensional and the Hamiltonian H3 has eight eigenvalues 
(see Appendix lB)l . shown in Fig.Q]as functions of the 
coupling constant A. Two of them are double degener- 
ated, while the remaining four arc not, apart from several 
isolated values of the parameter A. The ground state of 
the system for any finite value of A is non-degenerate, and 
in our notation it is the seventh state \ej). When the pa- 
rameter A is infinite, which corresponds to the zero value 
of the external magnetic field B, the Hamiltonian H3 has 
only one free parameter, and can be expressed as 

H 3 (X = 00) = -Ci (of <g> erf + erf <g> erf + of <g> erf) . 

The two lowest states (denoted as \gi) and \g2)) in the 
energy spectrum become degenerate in this case. These 
states read 

I.91) = 5(1000) + |110) + |011) + |101)), 
\gi) = ^(|iii) + |ioo) + |oio) + |ooi)). 

Here the first one is the limit of the state \ei) when A 
tends to infinity and the second one is the limit of the 
state I es). We know that any linear combination of the 
two vectors is an eigenvector with the same energy and, 
consequently, can be identified as a ground state. How- 
ever, there is one exceptional linear combination. If we 
demand the ground state of the system with the param- 
eter A = 00 to be the limit of the ground state when 



A = 00, then the appropriate choice for the ground state 
is |Si)- 

As we will see in Section TlV Al the point A = 1 turns 
out to be rather interesting. There is a particular eigen- 
state of the Hamiltonian which has quite interesting be- 
haviour of entanglement around A = 1. However, for N 
large it is rather difficult to identify this specific state 
among 2 N eigenstates of the Hamiltonian H^. Having 
calculated the spectrum, the state can be easily identified 
with the help of the level crossing at the point A = 1 . The 
special state with the remarkable properties is in general 
non-degenerate but at the point A = 1 becomes degener- 
ate 

E 8 (X = 1) = E 3A (X = 1) = , 

and crosses the degenerate levels E 3 and £4. What is 
important is the fact that this type of level crossing is 
independent of N (we might say universal) , there is the 
same type of level crossing for N being an arbitrary odd 
number. 



B. Entanglement properties 

Our main goal is to analyse the entanglement prop- 
erties of the model. Let us begin with the Ising chain 
of three qubits in the ground state and examine entan- 
glement as a function of the parameter A. We will use 
the three different measures of entanglement: The con- 
currence, the tangle and a measure of the intrinsic three- 
partite entanglement, as introduced in Section III Bl An 
important aspect is the comparison of specifically bipar- 
tite, and multi-partite entanglement. 

The bipartite entanglement between individual qubits, 
the entanglement between a qubit and the rest of a sys- 
tem, and an intrinsic three-partite entanglement for the 
ground state are shown in Fig. [21 They are quantified by 
concurrence, tangle, and the intrinsic three-partite en- 
tanglement of Eq. H2.6I) . respectively. Due to the fact 
that any non-degenerate state shares all symmetries of 
the corresponding Hamiltonian, the entanglement of the 
ground state between an arbitrary pair of qubits has the 
same dependence on the parameter A and it holds that 

C 12 (A) = C 13 (A) = C 23 (X) , (3.3) 

Moreover, the same holds for bipartite entanglement be- 
tween a given qubit and the rest of the system so that 

Gn(A) = C a2 (A) = Caa(A) , (3.4) 

where X denotes a system of two qubits with the qubit on 
the X-th position omitted and C X x 1S the entanglement 
shared between the qubit on the X-th position and the 
rest of the Ising chain (remaining two qubits). 

The solid line in the picture Fig. [21 shows bipartite 
entanglement between an arbitrary pair of qubits. For 
A = the concurrence, i.e. the pairwise entanglement is 
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2 4 ^6 8 10 

FIG. 2: The entanglement in the ground state of the Ising 
chain with three qubits JV = 3. In the figure we present 
three different types of the entanglement as a function of the 
parameter A: the bipartite entanglement between the first and 
the second qubit C12 (solid line); the bipartite entanglement 
between the first qubit and the remaining two qubits Cu 
(dotted line) and the intrinsic three-partite entanglement C123 
(dashed line). 



zero. As A = Ci/B it means that the case A = corre- 
sponds to the absence of the interaction, Ci = 0. Con- 
sequently, the ground state of the system is such that all 
spins are aligned along the same direction, the direction 
of the magnetic field, and are not entangled. When we 
turn on the interaction, the constant Cj is no longer zero 
and the spins become entangled. As we increase the value 
of the interaction constant Cj (or equivalently decrease 
the value of the magnetic field B so that the ratio Cj / ' B 
increases) the qubits become more and more entangled. 
This holds up to the value of A = 1 where the two qubit 
entanglement reaches its maximum. Further increase of 
the parameter Cj (or decrease of B) causes degradation 
of the entanglement and in the limit A goes to infinity 
the entanglement is zero. As we have chosen the ground 
state for A = 00 to be the limit of the ground state when 
A tends to infinity, the concurrence at the point A = 00 
is zero. It means that when the magnetic field is zero, 
the ground state of the system is such that all pair con- 
currencies are zero and there is no entanglement in any 
pair of qubits. 

The entanglement of a given qubit X and the rest of 
the system X expressed in terms of the tangle is, on con- 
trary, a non-decreasing function of A. At the point A = 
the entanglement is zero on account of the same reason 
as the entanglement between an arbitrary pair of qubits. 
When the parameter A is non-zero (i.e. the interaction 
constant Cj is non-zero) the qubits are entangled. That 
is, any chosen individual qubit is entangled with the rest 
of the system. The stronger the interaction (the larger 
the value of the parameter A), the stronger the qubits 
are entangled with the system. In the limit A — > 00 
(infinitely strong interaction) the qubits become maxi- 
mally entangled and the tangle, measuring the amount 
of entanglement between a given qubit and remaining two 
qubits, reaches its maximum value. For the case of our 
specific choice of the ground state for A = 00 the tangle 



is maximal and equals to unity. 

As we have already pointed out the reason why we have 
described the case of three qubits in such detail is that 
the entanglement behaves in the same manner for an ar- 
bitrarily large N. But the case of three qubits is special 
for a different reason too. In the case of just three qubits 
being in a pure state we are able to calculate the intrinsic 
three-partite entanglement using Eq. (|2.6I) . The dashed 
line in the Fig.[2]shows the dependence of the three-qubit 
entanglement on the parameter A. We can see that the 
dependence of the intrinsic three-partite entanglement on 
the parameter A is very similar to the dependence on 
the same parameter of the entanglement between a given 
qubit and remaining two qubits (Fig. [21 . It seems that 
for a strong interaction the three-partite entanglement is 
the largest contribution provided we express the entan- 
glement between a single qubit and remaining two qubits 
(rest of the system) as a sum of two and three-qubit en- 
tanglement (see Eq. (I2.6[l and comments therein). This 
result suggests the following physical picture: "When the 
system of interacting spin- 1/2 particles is in the ground 
state then the interaction causes entanglement of qubits 
such that each qubit is entangled with the rest of the sys- 
tem. For the system of N spin-1/2 particles the N -partite 
entanglement will be dominant when the interaction be- 
tween the particles is very strong compared to the mag- 
nitude of the magnetic field. " This conjecture, proven to 
be valid in the case of three qubits, will be further ex- 
amined in following sections where the general case of a 
chain with an arbitrary number of qubit will be analysed. 



C. Entanglement in Gibbs ensembles 

In this subsection we will continue to investigate the 
3-qubit Ising model. We will analyse the entanglement 
properties of thermal states of three qubits interacting 
according to the Ising Hamiltonian. The ground state 
of the system is probably the most important state and 
through the study of those states we acquire a lot of in- 
formation about the corresponding system itself. Beside 
being the states with the lowest energy we know that 
the ground states are associated with zero temperature 
and that they are related to fundamental properties of 
Hamiltonians. However, there are other states which are 
equally relevant for the physical description of the sys- 
tem. The temperature of a system we measure in our lab- 
oratories is always non-zero. Keeping in mind the third 
law of thermodynamics and the impossibility of reaching 
the absolute zero temperature we can conclude that in 
practice there is always a non-zero probability for find- 
ing the system under study in one of the excited states. 
Of course, the probability depends on the temperature 
but as far as the temperature is non-zero, no matter how 
big the gap in the energy between the ground state and 
the first excited state is, the probability is non-zero as 
well. Consequently, it is interesting to study entangle- 
ment in systems in thermal equilibrium, i.e. in their 
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FIG. 3: The bipartite entanglement between the 1st and the 
2nd qubit as a function of the temperature T and the param- 
eter A. The system of three qubits is in the thermal state 
Eq. ll33t . 



"natural" state, and investigate the dependence of the 
entanglement on temperature. 

The density operator corresponding to a thermal state 
of a quantum system at the temperature T is usually 
given by the relation 



P( T ) = ^2wi\ei)(ei\ 



(3.5) 



where |e<) is an energy eigenstate (eigenstate of the 
Hamiltonian H3), Wi are weights or probabilities defined 



= Ke 



-Ei/T 



(3.6) 



where we assume the Boltzmann constant to be equal 
to unity and we sum over all energy eigenstates. The 
constant K in Eq. I|3.6|l is a normalisation so the sum of 
probabilities Wi equals to unity 



= 1 



(3.7) 



In Fig.[5]we have plotted the entanglement between the 
first and second qubit in a three-qubit system at temper- 
ature T. Repeating the same line of arguments, taking 
into account the symmetry of the Gibbs state at the tem- 
perature T, we know that C\i = C13 = C23 and Fig. |3 
shows us the dependence of the entanglement on tem- 
perature for an arbitrary pair of qubits. For nearly zero 
values of the temperature the entanglement behaves in a 
similar way as in the case of the system in the ground 
state. Increasing the temperature the two qubits become 
less and less entangled. In the high temperature limit the 
entanglement is practically equal to zero. It has a very 
simple explanation. If the temperature is high enough 
all probabilities u>i are almost equal and the state of the 
system p is proportional to the identity operator [34J, i.e. 
it is the total mixture. Consequently, the state of an ar- 
bitrary pair of qubits is proportional to the identity as 
well and the two qubits are not entangled. 



In our case an increase of the temperature always 
causes degradation of the entanglement. Thus we may 
conclude that to maximalize the entanglement it is con- 
venient to keep the temperature as low as possible. It fol- 
lows that under certain conditions, one way of increasing 
the entanglement can be lowering of the temperature. At 
the end let us note that there are quantum models where 
an increase of the temperature can cause an increase of 
entanglement (see for instance Ref. Q). 



D. Quantum entanglement at A = 1 

Performing an analysis of the entanglement for the 
whole set of eigenstates of the Hamiltonian H3, we have 
found one particular eigenstate with rather peculiar be- 
haviour of the entanglement. Namely, the entanglement 
of this state as a function of the parameter A is non- 
analytic at the point A = 1. 

Let us remind the reader that at the end of Sec- 
tion [^^3 we have mentioned a level crossing. That is, 
at the point A = 1 there is an energy level crossing and 
one of the non-degenerate eigenstates becomes degener- 
ate. What is remarkable is the fact that the eigenstate 
with non-analytic behaviour of entanglement is the same 
state mentioned in Section IIII Al in connection with the 
level crossing. However, while we have discovered the 
state through our analysis of the entanglement of the 
eigenvectors for three qubits, in the general case of arbi- 
trary odd number of qubits we have followed the reverse 
path : we have identified the state by exploiting the level 
crossing at the point A = 1 [35^. 

In what follows we will denote the state exhibiting this 
very intriguing behaviour as the "X-state" (since it ex- 
hibits extreme entanglement around A — 1 - for details 
see Sec. IV. B). In our earlier notation, it is the eighth 
state \es). To remind the reader the state has the follow- 
ing form: 



\X) EE |« 8 ) = K 8 



1-E 8 -2X 
A 



|000) + |011) + |101) + |110) 

where Kg is a normalisation constant and Eg, is the en- 
ergy corresponding to the eigenstate |e 8 ) = \X). The 
X-state is a non-degenerate eigenstate of the Hamilto- 
nian (except for a finite number of values of A) and thus 
shares all symmetries of the Hamiltonian H3 in the same 
way as the ground state. That is C12 = C13 = C23 and 
Cyy = C 2 2 = C33- The bipartite entanglement between 
the first and second qubit C\i and between the first qubit 
and remaining two qubits C\\ are shown in Fig. 0] We 
from from the figure that the concurrence between two 
qubits in the system exhibit non-analytical behavior at 
A = l. 

Certainly, the reason behind this non-analyticity can- 
not be a phase transition. We know that the Ising model 
has a quantum phase transition at the point A = 1 but 
for that the chain must be infinite and the temperature 
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FIG. 4: The entanglement in the X-state of the Ising chain 
with three qubits TV = 3. In the figure we show two different 
types of entanglement as a function of the parameter A: the 
bipartite entanglement between the first and the second qubit 
C12 (solid line) and the bipartite entanglement between the 
first and remaining two qubits Cii (dotted line). 



must be zero |30j . That is, we can observe a phase transi- 
tion only if there is an infinite number of particles in the 
chain and the system must be in the ground state. From 
this point of view there cannot be a direct link between 
the found non-analyticity and the phase transition. 

The other question is the relation between the entan- 
glement and a change of symmetries in the system. The 
change of symmetries of a Hamiltonian can have a sig- 
nificant effect on the correlation properties of the eigen- 
states. In our case we know that the phase transition is 
accompanied with a symmetry change at the correspond- 
ing point. Similar change of symmetry is observed even 
in the case of a finite dimensional Ising chain at the point 
of the level crossing. This suggests that the change of the 
symmetry at the point A = 1 may in general be reflected 
in the behavior of the entanglement - the quantum part 
of the correlations. 

The X-state is interesting not only on account of the 
non-analyticity but mainly for the fact that it exhibits 
remarkable quantum correlations. At the point A = 1 
the entanglement in an arbitrary pair of qubits is zero 
but the entanglement between a given qubit and the re- 
maining two qubits is maximal (c.f. Fig. 0J. If we cal- 
culate the entanglement length introduced in Ref. pj, it 
is zero. But we know that each qubit is maximally en- 
tangled with the rest of the system so there is a sort 
of a "long range" entanglement [3(j. In other words, as 
each qubit is maximally entangled with the rest of the 
system and the pairwise entanglement is zero, we have 
an intrinsic multi-partite entanglement. Moreover, since 
the system consists of only three qubits the only possi- 
ble multi-partite entanglement is a three-partite entan- 
glement C123 and at the point A = 1 the three-partite 
entanglement reaches its maximal possible value. 

Tu sum up we can conclude that in the case of the finite 
dimensional Ising model, there is an energy eigenstate, 
the X-state, for which entanglement exhibits a rather 
special behaviour at the point, where the infinite Ising 
chain has a phase transition. 



IV. GENERAL CASE OF TV SPINS 

So far we have considered the particular case of three 
qubits. Despite its simplicity the case of three qubits 
shares many features of the general case of a chain with 
an arbitrary number of qubits. This has helped us to 
formulate basic theorems and to identify states which are 
particularly interesting with respect to the entanglement. 

Let us consider a chain with TV qubits where TV is ar- 
bitrarily large. The Hilbert space of TV qubits is 2 N - 
dimensional and the corresponding Hamiltonian Hn has 
2 N eigenvectors and eigenvalues. Despite the very pos- 
sibility to calculate any eigenvector or eigenvalue (recall 
that the model under study is exactly solvable) it is not 
feasible to perform the calculation for all eigenvectors 
(eigenvalues) and to analyse them afterwards. There- 
fore, we have used the results of the previous section and 
beside the ground state as an important state, we have 
analysed the X-state. Of course, prior to that we have to 
find or identify the X-state among 2 N eigenstates of the 
Hamiltonian H^. At this point we can take the advan- 
tage of our detailed knowledge on the spectrum we have 
at hand, as the level crossing at the point A = 1 studied 
in the previous section is crucial in identifying the sought 
state. 

Our main goal is to analyse the entanglement proper- 
ties of the states under consideration. In addition to the 
dependence on the number of qubits TV in the chain we 
will also consider the dependence on the "distance" of 
qubits. For TV > 3 there are more possibilities how to 
create pairs of qubits and beside the nearest neighbours, 
a pair can be created from the next nearest neighbours 
and etc. Since the Ising model is not permutationally 
invariant unlike in the special case studied so far, we can 
expect that the entanglement will vary with the distance 
between qubits. 



A. The ground state 

The ground state of the system for different values of TV 
can be calculated using the formalism developed in Ap- 
pendix lAl and Appendix IU1 Due to the complicated form 
of the state itself, we do not quote the explicit expression, 
it may be found in the above mentioned Appendices. In 
the following, we describe the entanglement properties of 
the state instead. 

The entanglement shared between pairs of nearest 
neighbour qubits in terms of the concurrence is plotted 
in Fig.[S]for different values of TV. The shape of different 
curves, corresponding to different number of qubits in the 
chain, is very similar. At the point A = the values of all 
curves are zero, increasing the parameter A the entangle- 
ment (measured in the concurrence) increases and around 
the point A ~ 1 reaches a maximum. This maximum de- 
pends on the number of qubits TV, but with increasing TV 
the concurrence converges to a specific value that even 
for TV > 5 is almost constant. Increasing the parame- 
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FIG. 5: The entanglement between the nearest neighbours as 
a function of the parameter A and the number of qubits in the 
chain N. The system is in the ground state and the number 
of qubits in the chain is N = 3, 5, 7. 



tcr A further, the entanglement decreases and finally, in 
the limit A — > oo the entanglement tends to zero. For 
A — ► oo, the ground state is degenerate. Similarly, as in 
the case N = 3, we may choose a particular ground state 
the A — ► oo limit of which becomes the ground state of 
the Hamiltonian for A — > oo (c.f Appendix [CJ 



JV/A=oo 



= no = K 



N 



{i,j,...}e 



{'•./••••I ••• 



(4.1) 



where Kjy is a normalisation constant, {i,j, . . .} denote 
positions of the qubits that are up and . . .}e means 
summation over all states of the standard basis with an 
even number of qubits up. (We use the term "up" for 
a qubit if it is in the state |1) and down if it is in the 
state |0)). From the construction of the state \iPn)\=oo, 
it follows that the entanglement between arbitrary two 
qubits is zero while the entanglement between a given 
qubit and all remaining qubits is maximal. 

Proof: The state Eq. (|4.1|l can be rewritten into a sim- 
pler form using the following recurrence relation: 



|V ) JV)A=00 = -J= |0)|^at_i) A =oo + |l>|^JV-l)A=oo(|4,2) 



where \ipN-i)\=oc has the same form as |^jv_ i)a=oo , but 
instead of summing over all states with an even number 
of qubits up we sum over all states with an odd number 
of qubits up. With the help of the Eq. (|4.2() it is easy 
to prove the above statements concerning entanglement. 
Let i and j denote two arbitrary but mutually different 
(i 7^ j) positions of a priori chosen qubits in the chain. 
Using the relation Eq. (|4.2II we rewrite the state vector 
Eq. 14.111 as follows: 



\<Pn)x=oc = \ [ ( |00)y + |ll)y )|^Ar_ 2 ) A= c 
+ (|01)tf + |10)y)|^_ 2 )A=« 



(4.3) 



The reduced density operator p^ , of the two qubits at the 
positions i and j is calculated as a trace over all remaining 



qubits of the density operator p — |V , jv)a=oc('0aH of the 
whole system. The state p^ expressed as a matrix in the 
basis {|00), |01), |10), |11)} reads 



(\ w 



Pij 



1 



(4.4) 



The spin-flipped density matrix in this case reads pij — p 
so that Pij Pij = \pi j- The matrix pijp%j has two eigen- 
values that are equal and according to Eq. 12.41 the en- 
tanglement shared by the two qubits in the z-th and j-th 
positions is zero. 

The density operator of a single qubit at the z-th po- 
sition pi can be calculated similarly, 



Pi 



(4.5) 



and it corresponds to a maximally mixed state. Since 
the whole system is in a pure state, the qubit in the 
z-th position is maximally entangled with the rest of the 
system. 

Let us note, that the state given by Eq. (|4.1|) is a GHZ 
state in the basis built by direct products of eigenvectors 
of o~ x of each qubit. Thus the state can be transformed 
via local unitary transformations into the standard form 
of the GHZ state 



\GHZ) 



1 



V2 



(|00...0) + |11...» 



(4.6) 



This observation provides us with a very clear explana- 
tion of the above mentioned entanglement properties. 
Moreover, it thoroughly confirms the proposal conjec- 
tured at the end of Subsection IIII Bl For A being large 
but not infinite or, equivalently, for a large coupling con- 
stant Cj compared to the absolute value of the magnetic 
field B, the ground state of the system exhibits properties 
close to the GHZ state. Taking the parameter A larger 
and larger the state is closer and closer to the GHZ state 
and for a sufficiently large A we can consider the ground 
state of the system to be the GHZ state even for A be- 
ing large but finite. When the state of the system is the 
GHZ state, the entanglement between any pair of qubits 
is zero because the reduced density operator describes a 
separable state [c.f. Eq. (|4.4() ]. Further, a reduced den- 
sity operator of a subsystem consisting of n < N qubits 
is also separable, as one would expect for a GHZ state. It 
follows that if we consider an arbitrary subsystem there 
is no entanglement: choosing any set of n < N qubits, 
the state of the chosen set is separable. Consequently, 
the state under consideration exhibits only intrinsic TV- 
partite entanglement. Recalling the conjecture from Sub- 
section IIII Bl we can now confirm the result to be valid 
for a general case of the iV-partite chain as well. 

Finally, let us note that the ground state (or equiv- 
alently the GHZ state) has a long range entanglement. 
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FIG. 6: The bipartite entanglement shared between qubits, 
expressed in terms of the concurrence dj where i and j label 
positions of two qubits, as a function of the parameter A and 
the distance between the qubits. The system (chain) of seven 
qubits is in the ground state. 



The TV-partite entanglement is certainly of a long range 
since it concerns all qubits in the chain, though, for in- 
stance, the entanglement length defined in Ref. [lj is zero. 

We have also studied how the entanglement depends on 
the distance between an arbitrary pair of qubits qubits, 
see Fig. The farther the two qubits are they are less 
entangled. It means that with increasing the distance 
(specified by positions of qubits, i.e. the distance between 
qubits ji and j'2 is represented by the difference \ji — 
j2 1 ) between qubits the entanglement converges to zero. 
Besides, the peak, or the point where the entanglement 
is maximal, is shifted to larger A's (see Fig. |SJ). 



B. General X-state 

In the case of the chain composed of three qubits we 
have found a particular eigenstate of the Hamiltonian H3, 
the X-state, with a non-analytic behaviour of entangle- 
ment. The question naturally arises, whether there exists 
such a state in the case of more than three qubits. As 
the Hamiltonian Hm has 2^ eigenvectors and eigenval- 
ues it is impossible to analyse the whole spectrum. How- 
ever, we know that in the three qubit case the X-state 
is interesting not only on the account of entanglement 
but also because of energy-level crossing. It is a non- 
degenerate eigenvector of H3 apart from a single point 
A = 1 where there is a level crossing. Our knowledge of 
the level crossing can be successfully exploited in identi- 
fying the X-state for an arbitrary N. We have found that 
for N odd there is a level crossing at the point A = 1 
and one of non-degenerate eigenvectors becomes degen- 
erate. Let us note that in the case N = 3 the X-state 
corresponds to the non-degenerate eigenvector while the 
other energy level is degenerate and corresponds to two 
vectors. The situation we have now is similar. There is 
a level crossing of two energy levels at the point A = 1. 
One of them is a non-degenerate energy level; in what 
follow we will call the state corresponding to that level 
as the X-state, and the other energy level is degenerate 



and there are 2™ (where N = 2n + 1) eigenvectors cor- 
responding to that level (for proofs and more details see 
Appendix [DJ . 

Having successfully identified the X-state, we can anal- 
yse its entanglement properties. In contrast with the case 
of N — 3 qubits, for N > 3 bipartite entanglement (the 
concurrence) as a function of the parameter A is analytic. 
In order to see this we need to know the form of the X- 
state. It is derived in Appendix ^ we quote only the 
result for A = 1 here (see Appendix IU1 and also Ref.jU) 



I*) 



A=l 



E 

{lj, ■■■}<! 



\{i,j,k,...})(-l)^>i 



(4.7) 



where {i,j,...} e denote a sum over all sets of indices 
with an even number of indices in each set, the letters 
i, j, ... in a single set denote positions of qubits in the 
chain that are up and •••}) is the corresponding 

state vector, d(i,j) is a distance between qubits on the 
i-th and j-th positions defined below, while the sum over 
i < j means that we sum over all pairs of qubits counting 
only once the cases with switched positions of the qubits. 
The distance d(i,j) of the two qubits is defined as the 
shortest path on the ring that brings us from the qubit 
on the i-th positions to the qubit on the j-th position. 
In Appendix [D] we present a complete proof that the 
state Eq. (|4.7() is an eigenstate of the Hamiltonian Hn for 
A = 1. The Appendix iDl also contains several additional 
proofs and more details on the A-state. 

The reduced density operator of two qubits on the i- 
th and j-th positions has been calculated from the state 
Eq. H4.7fl by tracing over degrees of freedom of the re- 
maining qubits. The density operator expressed in the 
basis {|00), |01), |10), |11)} reads 



pij = 



(\ \ 

\ 






\ 



(4.8) 



\o \) 



In the same way the density operator of a single qubit 
in the i-th position j38| expressed in the one-qubit basis 
{|0),|l)}is 



Pi 



* ? 
h 



(4.9) 



With the help of the density matrix Eq. i|4.8|) we have 
calculated bipartite entanglement between qubits on the 
i-th and j-th positions while the density matrix Eq. Ij4.9|l 
has been used for the calculation of entanglement 
between a qubit on the i-th position and the rest of 
the system (all remaining qubits). Since the density 
matrix Eq. (|4.8Jl is a complete mixture, there is no 
bipartite entanglement between any two qubits. On the 
other hand as Eq. (|4.9() describes a complete mixture 
and the whole system is in a pure state, entanglement 
between a given qubit and remaining qubits is obviously 
maximal: the tangle is equal to one. As the eigenvalues 
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of the density matrices pij and pi are continuous 
functions of the parameter A it is easy to check that 
both the entanglement shared between qubits and the 
entanglement between a given qubit and the rest of 
the system are continuous functions of the parameter 
A around the point A = 1. Consequently, the point 
A = 1 is not a point of non-analytical behaviour of 
entanglement anymore. Moreover, the state Eq. (|4.7|) 
is not equivalent to a GHZ state. In other words there 
does not exist a local unitary transformation which 
would transform the state Eq. 1)4. 7fl into the GHZ 
state in Eq. I|4.6[l . However, we have found that the 
X-state Eq. (|4.7() has the following remarkable property: 

Theorem 1 Let N = 2n+l denote the total number of 
qubits forming the Ising chain where n is an integer and 
let the system of N qubits be in the X-state Eq. \4- 7| ). 

The density matrix of any sequence of n neighboring 
qubits is 

P=^*n, (4-10) 

where l n is the identity operator acting in the 2 n - 
dimensional Hilbert space of n qubits. 

An important consequence of the theorem is the fact that 
if the system is in the X-state then any set of neighboring 
qubits is perfectly (maximally) entangled with the rest of 
the system. Consider an arbitrary set of n neighboring 
qubits. The reduced density operator of the n qubits is 
Eq. (|4.10|l . Using the entropy Eq. 1)2. 7JI as a measure of 
bipartite entanglement for pure states 

s = -E^ lo 4 =,llog2 ' (4 - n) 

2" 

we can see that the set of n qubits is maximally entan- 
gled with the remaining qubits. That is, if a system of 
N = 2n + 1 qubits is in the X-state and we choose n 
neighboring qubits then according to Eq. (|4.11|l we know 
that the n qubits are perfectly entangled with the re- 
maining n + 1 qubits. Moreover, if we choose a subset 
of say k qubits from the set of n neighboring qubits then 
the state of the k qubits is 

P=^lfc> (4-12) 

so the k qubits are perfectly entangled with the rest of 
the system (all remaining N — k qubits). To sum up the 
X-state is a highly entangled state and consequently it is 
a good candidate for a quantum communication between 
many parties. A rather simple protocol that can serve 
as an example of its applications described in the next 
section. 



C. Controlling distribution of entanglement in the 
X-state 

We will present a simple example how the X-state can 
be exploited for a communication or a secret key distribu- 
tion in a situation when bipartite entanglement between 
many parties is needed. The X-state with its remarkable 
properties can be considered to be a very good resource 
of communication as any set of n neighboring qubits is 
maximally entangled with the rest of the system. 

First, imagine that n neighboring qubits belong to Al- 
ice and the rest (n+ 1 qubits) belongs to Bob. Moreover, 
let us assume that Alice and Bob want to exploit the en- 
tanglement of the X-state for their protocol. But unfor- 
tunately, their protocol is designed for qubits, that is to 
say it uses pairs of maximally entangled qubits. We have 
shown that the density operator of any pair of qubits is 
proportional to the identity and thus the two qubits can- 
not be entangled. It means that Alice and Bob cannot 
take any two qubits and use them for their protocol. But 
we know that the n neighboring qubits which belongs to 
Alice are maximally entangled with Bob's qubits as the 
entropy in Eq. (|4.11|l equals to n log 2. Such amount 
of entanglement corresponds to n pairs of maximally en- 
tangled qubits. Therefore one may ask whether they are 
able to create n pairs of maximally entangled qubits only 
by performing local (though multi-qubit) unitary trans- 
formations and U^ l+1 ' ) on their respective qubits. 
The answer is positive. 

Consider a state of 2n + 1 qubits with n pairs of maxi- 
mally entangled qubits and let Bob's last qubit be in the 
state |0) (We know that unlike Alice, Bob has got n + 1 
qubits.): 




Now let us reorder the qubits in such a way that the 
first n qubits belong to Alice and the remaining n + 1 
qubits belong to Bob. We need to do that because Alice 
possesses n neighboring qubits. 

|E> = J2 |{U,--.}>^i|{U,-.J>®|0>(4.13) 

where a set of indices {i, j, . . .} denotes positions of qubits 
up in the standard basis vector and we sum over all pos- 
sible sets of indices. Now we want to find a local unitary 
transformation U = U^* 1 ' such that the state 
|H) transforms into X-state. It follows from Eq. i|D13|) 
that it is enough to consider the unitary transformation 
of the form U = 1^ <g> Ug +1 ' where 

U& +1) :K iJ ...})o-^7jlKj,---}>®|0), 

and the states \otuj..\)o are defined in the Appendix iDl 
After Bob has performed the unitary operation Alice and 
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Bob share n pairs of maximally entangled qubits and they 
can begin with their original protocol. This simple exam- 
ple illuminates the remarkable properties of the X-state 
and its convenient form since only Bob has to perform 
the local unitary operation. 

The situation becomes even more interesting if we re- 
place Alice with n parties {Ai, . . . A n }. Now, Bob com- 
municates with n different parties. By performing a local 
operation he can decide which of his qubits is entangled 
with a given partner Aj . Let us stress that this is only a 
simple example and more sophisticated protocols are the 
topic of current research. 



V. SUMMARY 

We have performed detailed analytical calculations 
concerning stationary states of a finite-size Ising chain 
with cyclic boundary conditions, and their entanglement 
properties. We have put a special emphasis to a kind of 
description of multi-partite entanglement. 

The primary motivation of our investigation has been 
an attempt to illuminate the Ising model using tools of 
quantum information theory. In addition we were study- 
ing a possibility whether some properties of the Ising 
model can be used as a resource of quantum informa- 
tion processing / communication. For this purpose, one 
of the criteria that should be met is that of the possi- 
bility of preparing the system in a suitable initial state. 
As physical systems tend to occupy their ground states, 
it is advantageous if the ground state is a suitable initial 
state for some purposes We have shown that by adjusting 
the external magnetic field the ground state of the model 
considered is the well known GHZ state used in several 
quantum information processing schemes. Consequently, 
with the ground state of the system well known and hav- 
ing particularly nice properties makes the Ising model a 
good candidate as a resource for a quantum information 
processing. 

This result also demonstrates the usefulness of the ap- 
proach to finding an entangled state with pre-defined 
multi-partite entanglement properties by finding the 
ground state of a suitably chosen Hamiltonian. The 
ground state of the Ising model for certain values of the 
parameter A is a very specific state - the GHZ state. The 
GHZ state has the property that the entanglement be- 
tween any set of n qubits where n < N is zero while 
the iV-partite entanglement peaks reaching the maximum 
possible value in the limit A — > oo. It means that in 
the limit A — > oo the ground state of the Ising model 
maximises the iV-partite entanglement and the state of 
N qubits with maximum A-partite entanglement can be 
found as a ground state of the Ising model with N qubits 
in the chain. 

Regarding entanglement properties, not only the 
ground state of the Ising Hamiltonian is found to be in- 
teresting. For instance, we have discovered a very spe- 
cific eigenstate of the Hamiltonian, the X-state. The X- 



state is strongly (extremely) entangled as every set of 
k < n neighboring qubits (where the total number of 
qubits is N = 2n + 1) is maximally entangled with the 
remaining qubits. An important message of our results 
is, that multi-partite entanglement plays a crucial role 
in the understanding of exactly solvable models of quan- 
tum statistics. We have also presented a simple example 
to demonstrate the usefulness of such state in quantum 
communication. 

The X-state is identified via a certain crossing of en- 
ergy levels at A = 1, where a phase transition occurs in 
the thermodynamic (infinite-qubit) limit. Consequently, 
there might be some connections between a functional de- 
pendence of the entanglement as a function of A and the 
point of a phase transition. One important lesson one can 
learn from our investigation is that higher energy eigen- 
states of the Ising Hamiltonian might carry non-trivial 
information about quantum correlations of the system 
under consideration. 

This was work supported in part by the Euro- 
pean Union projects QGATES, QUPRODIS, and CON- 
QUEST and by the Slovak Academy of Sciences. We 
would like to thank Matyas Koniorczyk for his helpful 
comments. VB thanks Michael Nielsen for stimulating 
discussions and hospitality during his visit to the Queens- 
land University in 2000 when this work has been done. 



APPENDIX A: EXACT SOLUTION OF THE ID 
ISING CHAIN 

The Ising model as one of the simplest exactly solvable 
models has been studied extensively in the literature and 
there is a chapter on the Ising model in almost every text- 
book on solid state physics with exactly solvable models. 
The derivation here will mainly, up to a few minor devi- 
ations, follow Ref. ,32-j . The reason why we summarise 
the derivation here is the fact that we present a complete 
solution together with the exact form of the eigenvectors 
usually omitted in the literature and also to keep our 
discussion self-content. 

Consider the Hamiltonian Eq. (|2.3|) with the cyclic con- 
ditions Eq. (|2.2() . Our task is to find the eigenvalues and 
eigenvectors of this Hamiltonian. As a first step we per- 
form several transformations. This is a bit technical part 
however it needs to be included in the derivation to keep 
it transparent. 



1. Transition from variables a a to a 

Firstly we introduce new variables where 
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The Hamiltonian Eq. I|2.3(l after the first transformation 
reads 

N 

H = + < J ™)( cr n+l +°n+l) 

n=l 

W 1 

+2^( ( t>---1) . 



A" 



2. Jordan- Weyl transformation 

Secondly, we introduce fermionic variables <y n and Cn 
such that 

n-l 

c„ = exp(i7r^£7-to-7)c7-- , 

3=1 
n-l 

4 = ex P(* 7r Xl Cr / ,T 3 7 ) Cr « ■ 

3=1 

The variables cjj and c„ satisfy the anti-commutation re- 
lations. The form of the Hamiltonian after the second 
transformation is following 

N-l 

H = -A ((4 - Cn)(4n+1 + C ™+l) 

n=l 

-(-l)*(4_ Cjv )( c t +c ^ 

N ( 1 

n=l V J 

where N — X^=i 4 C « can De interpreted as the operator 
of number of fermions. The operator (— 1) commutes 
with the Hamiltonian H . For that reason it is possible 
to choose common eigenvectors of the Hamiltonian and 
the operator (— 1)^. The eigenvalues of the operator N 
are even or odd. If an eigenvector of the operator N cor- 
responds to an even eigenvalue then the operators cn+i 
and cjy +1 satisfy the condition 



cn+i 



anticyclic cond. 



-Ci . 



On the other hand if an eigenvector of the operator N 
corresponds to an odd eigenvalue, the operators cjv+i and 



cj v+1 are defined as 



J 
-l 

XC N+ i = Ci 



cyclic cond. 



The advantage of the previous choice is that the Hamil- 
tonian has the same form in both cases 



N 

H = "A ^ ((4 - C ")(4+l +Cn+1. 



+ 2 £ 4 



n=l 



3. Momentum representation 

The last transformation is a transition to the momen- 
tum representation and specified by the variables rf q and 
r\ q . The fermionic operators c' n and c„ and the new op- 
erators and r\ q are related as 



= —y 

.Sat ^ 



N 

-y 

q 



e iqm -n q 



e- iqm r]\ 



where q depends on the boundary conditions. For cyclic 
boundary conditions (TV is odd) we have 



e iqN = 1, 



so that 



N 



21 l = 0,l,2...(N -1) 



For anti-cyclic boundary conditions (iV is even) 
e iqN = -1, 

we have 

q = ^(21 + 1) 1 = 0,1,2. ..(N-l). 

The Hamiltonian H has again the same form in both 
cases and is given by 

-Xi sin q{vlvL- q + r] q r}2^- q ) 
+2 Hv q - \) 

It is important to note that the operator TV can easily be 
rewritten using the new variables rj q and rf q 



N = E c « c " = • 



(Al) 



It follows that the number of fermions described with 
c m and is the same as the number of other fermions 
described with the operators r\ q and rf q . 



4. Eigenvalues and eigenvectors of H q 

In what follows we calculate eigenvalues and eigenvec- 
tors of the Hamiltonian H. Firstly, it is convenient to 
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reorder the contributions to the sum Eq. (|A1|) . Specifi- 
cally, choose one q, and let us define H q (we note, that it 
is not possible to perform this step when q = or q = n. 
In that case see the end of this paragraph.) 

H q = -Acosg(2^r? g + 24 w _g%7r-g) 
-Xi sinq(2r]lj]l 7r _ q + Iri^-q) 

+ 2 {V% +vL-q 7 l27T- q - 1) • 

The Hamiltonian H is a sum of H a 



H 



At this moment we postpone the derivation of the eigen- 
vectors and eigenstates of the Hamiltonian H to the next 
paragraph and instead we discuss the cases when q = 
(or q = 7r) 39] . The form of operator Ho ( H n ) is 

H Q = -2 A ril ?/o + 2 eta Q - 1 q = 

( H„ = +2X^ + 2^-1 q = <K ). 

In this case we define the basis as follows 



|*,=o) = 77510 > 
|*o) = |0> 



[ !<!>.>= r?t|0 > ] 
[ |*o> = |0> ] . 



where we sum over half of the q's since H q = H2Tx- q . The 
point in rewriting the Hamiltonian as a sum of H q is that 
it is easy to calculate the eigenstates and eigenvalues of 
H q . For each H q we define the following basis 



= ri\o> 



l* 2 ) = nL- q \Q > , 

l$ 3 > = |o>, 

1*4) = vL- q vt\°>- 

In this particular basis the Hamiltonian H q is a 4 x 4 
matrix 



/-2Acosg 

-2Acos<? 

0-2 -2Aising 

\ 2Xisinq 2 — 4A cos q 



and the four eigenvalues of H q are 

a\ = — 2Acosg , 

a-2 = — 2Acosg , 

a s = 2(-A cos q + \J (A - l) 2 + 2A(1 - cos qj) 

a 4 = 2(-A cos q - y/{X - l) 2 + 2A(1 - cos qj) 



The eigenvectors corresponding to the eigenvalues <ii, 02, 
03 and a 4 in the defined basis reads 

|02> = 1*2) , 

|a3> = rfal^s) +e 3 |$ 4 ) , 
|a 4 ) = d 4 |<J> 3 ) +e 4 |$ 4 ) , 



where 



2 A i sin g 

1 



The operator f/o (-f^ ) m this particular basis is a 2 x 2 
matrix 



-1 
-2A+1 



The matrix has two eigenvalues 

<2l = — 1 

a 2 = -2A+ 1 



-1 
2A+ 1 



Oi = -1 ] , 
a 2 = 2A + 1 ] , 



with the corresponding eigenvectors given by 



= l*o) 
|oa> = |$,=o) 



[ M - l*o) ] 

[ |02) = I*.) ] 



5. Eigenvalues and eigenvectors of H: Example 

According to the previous paragraph the Hamiltonian 
H can be rewritten as a sum of the operators H q . Since 
each operator H q is a sum of terms with an even number 
of fermionic operators and two different H^s contain dif- 
ferent fermionic operators we conclude that the operators 
H q commute. Essentially it means that they have com- 
mon eigenvectors. As an example let us consider the case 
N = 5 and let (— 1) N . Then one of the energy eigenstates 
is 

\ e n) =^0^/5^/51°) 
with the corresponding eigenvalue 

2tT 47T 

E H = -1- 2Acos( — ) - 2Acos( — ). 

5 5 

We have to keep in mind here that not all the com- 
binations are the eigenstates of H. For example if 
(-1)* = -1 then the state 



\<P) = r L/B r L/B\°) > 



(A2) 



-, . 4A 2 (sing) 2 
1 + (a 3 +2) 2 



is not an eigenstate of H even though this state is an 
eigenvector of Hq, #2-7175 and H 47T / 5 . The state Eq. (|A2|) 
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has even number of fermions and consequently for TV even 
the Hamiltonian is not 



parameter A the Hamiltonian iJ/v has the following form: 



but 



H — H o + H 2rr /5 + H ijr / 5 , 



Recall that g's depend on the boundary conditions which 
in turn depend on the number of fermions. 



APPENDIX B: THREE SPINS: EIGENVALUES 
AND EIGENVECTORS OF THE HAMILTONIAN 

H 3 . 

For completeness we review the complete spectrum of 
the Hamiltonian H3. The Hubert space corresponding to 
the system of three qubits is eight dimensional and the 
Hamiltonian H3 has eight eigenvalues, 



-Ei, 2 

-^3,4 

E 6 
E 7 

Ee 



= A + l , 
= A-l, 

= 1-A-2V1 + A + A 2 



= 1-A + 2V1 + A + A 2 , 
A - 2^1-A + A 2 



-1 



-I-A + 2VI-A + A 2 



The eigenvalues E\ = E2 and E3 = E4 are degenerate for 
all values of parameter A while the other four are (apart 
from a finite number of points) not. The eigenvectors 
corresponding to the eigenvalues read 



e 2 



P3 



('i 



('a 



e-7 



1 

71 
1 

1 

71 
1 



[|110>-|011>] 
[|101> - |011>] 
[|100> - looi>] 



= _ [|010) - |001)] 



= K« 



K 7 



= K* 



E 5 


+ 1- 




E 5 


-3- 




Eq 


+ 1- 




Eq 


-3- 




1 - 


E 7 - 


2A 




A 




1 - 


Ea- 


2A 



111) 



111) 



1 000) 
1 000) 



|ooi> -i 

|001) 4 
f |011) 
f |011) 



|010) + 1 100) 
|010) + |100) 
f |101) + |110) 
f |101) + |110) 



N 



'i+i ■ 



The ground state of the Hamiltonian is degenerate, the 
energy of the ground state is E g — —NCi and two energy 
states related to the energy E g = —NCi are 



\9i) 

\g 2 ) 



|oo. 

= 111 



(Cl) 



where |0) and |1) are eigenstates of <j x . Of course not 
only these two states are eigenstates of the Hamiltonian. 
Any linear combinations of these states is also an eigen- 
state. However, we are searching for eigenstates that are 
the limit of the eigenstates of the Hamiltonian l|2.3|l when 
A tends to infinity. At this point we use the knowledge 
of the exact solution (see Appendix lA^I . and more specifi- 
cally the fact that any eigenstate of the Hamiltonian is a 
linear combination of standard basis vectors with cither 
even or odd number of qubits up. (Standard basis cor- 
responds to state vectors that are eigenstates of all erf. 
The fact that a state is a linear combination of standard 
basis vectors with either odd or even number of qubits up 
is equivalent to the condition that the state is an eigen- 
state of the operator (— 1) N .) Therefore we select two 
particular linear combinations, such that they are sums 
of state vectors of standard basis with either even or odd 
number of qubits up. The linear combinations of the two 
vectors in Eq. I|C1|I which satisfy the condition are 

1 

vT 
1 

V2 V 

where the first one is a linear combination of standard 
basis vectors with an even number of qubits up and the 
second one is a linear combination of standard basis vec- 
tors with an odd number of qubits up 

\gi)=K N Yj •••». 

{i,j,-}e 

\92)=K N IK Or ..•}>, 

{i,j,...}o 

and the constant Kn is a normalisation constant. It is 
easy to show that these are the only two possible linear 
combinations that satisfy the condition, and on top of 
that it is clear from the construction that each of the 
states |<7i) and I32) is a GHZ state. 



APPENDIX D: THE X-STATE 



\gi) 



\92) 



= (|00. 
= (|00. 



■> + |lT. 



Ill 



APPENDIX C: GROUND STATE FOR A = 00 

The case A = 00 corresponds to the physical situation 
with zero magnetic field B. For the given value of the 



The expression for the X-state reads 

\X)= |{M,^...})(-l) E *>^ (iJ) , (Dl) 

{i,3,...}e 
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FIG. 7: The ring of 9 qubits. The arrows denote two possible 
paths from the first qubit to the seventh qubit. The shorter 
path is the distance between the two qubits d(l, 7). 



where {i, j, . . .} e is a set of indices with an even number 
of indices in the set so that we sum over all sets of indices 
with an even number of indices, the letters i, j, . . . denote 
positions of the qubits in the chain that are up, d(i,j) is 
the distance of the qubits on the i-th and j-th positions 
defined below and the sum over i < j means that we sum 
over all pairs of qubits up counting only once the cases 
with switched positions of the qubits. Let us note that 
the state in Eq. (|D1I) is not normalised to unity. 

The distance d(i,j) of the two qubits is defined as the 
shortest path on the ring that brings us from the qubit 
on the i-th positions to the qubit on the j-th position. As 
the qubits form a circle, there are always two paths we 
can go without going backward and we can always choose 
the shortest one. To make clear what the distance defined 
above is, let us have a look at a simple example. Let N 
be 9 so that overall number of qubits in the ring is nine 
as in Fig. [7| Further let i = 1 so that it denotes the 
first position and j — 7 so that it denotes the seventh 
position. Then the shortest path is going from first to 
ninth position as the two are neighboring positions then 
from the ninth to the eighth and finally from the eighth to 
the seventh position. Consequently, the distance d(l, 7) 
in this particular case is 3. 



1. Proof that the X-state is an eigenstate of the 
Hamiltonian with zero energy 

In what follows we show that the state is the eigenstate 
of the Hamiltonian with zero energy or equivalently, that 
the following relation holds 



If we divide the Hamiltonian into the free Hamiltonian 
and the interaction Hamiltonian H = Hi + Hq, then the 
last equation can be rewritten as 

Hj\X) = -H \X) . (D2) 

The task now is to show that the two vectors - one on the 
left and the other on the right side of the last equation are 
equal. As we know, the equality of two vectors follows 
from the equality of their components in any complete 
basis. Actually what we will prove is the equality of 
the components of the two vectors in the standard basis 
(computational basis) i.e. 

{H I \X)) i = -(H \X)) i . (D3) 

Take one vector of the standard basis that is included 
in the sum given by Eq. (|D1J| and denote it as \v). We 
show that the v-th components obey Eq. (|D3(I 40] 

a. All vectors of the standard basis are eigenvectors 
of H . If we denote K to be the number of qubits up in 
the vector \v) then 

H \v) = (2K - N)\v) , 

and the v-th component of Hq\X) is 

(H \X)) V = (2K-N).s, (D4) 

where s is the phase factor of the vector \v) in the sum 
in Eq. (jDljl . 

b. Now it remains to show that the same holds for 
Hj except for the sign that must be opposite. The Hamil- 
tonian Hi is a sum of many elements C; where 

If we want to count the v-th component of Hi\X) we 
need to know the individual contributions from each term 
Ci\X). What is the action of the operator C;? It flips 
two neighboring spins on the j-th and (i+ l)-th positions. 
Let us assume that there are K spins up in the vector. 
If the two spins on the i-th and (i + l)-th positions are 
up then the operator C; flips them down and there are 
K — 2 spins up in the vector Ci\v). Similarly if the two 
spins on the i-th and (i + l)-th positions are down then 
the operator flips them so that they are up and con- 
sequently, there are K + 2 spins up in the vector Ci|w). 
Otherwise if one spin is up and the other is down then 
the number of spins up in the vector Ci\v) equals K. At 
this point it is obvious that neither Hq nor Hi can pro- 
duce a contribution (vector of the standard basis) with 
an odd number of qubits up since in Eq. 1|D1|) we sum 
over all sets of indices with an even number elements in 
each set. Now we use a little trick, namely 

C i = {v't a i+i) 2 = 1 1 

so that 



H\X) = 
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in order to answer the question which vectors from the 
sum in Eq. 1JJ1D give contributions to the v-th element 
considering only one Cj. Using the last relation the only 
possible one is 

-Sid\v), (D5) 

where &j is the coefficient of the state —Ci\v) in the sum 
in Eq. (|E> If) (We have introduced the minus sign in the 
last equation because of the sign of the operator C, in 
Hi.). Since the coefficients of the vectors in the sum 
Eq. (|D1|I are plus or minus one there are exactly N con- 
tributions to the v-th element, as we have N operators 
Cj, and all of them are plus or minus one. Our task is 
to find out the sign of each individual contribution and 
count them. First we divide the contributions into three 
subsets. Let us denote by Mo the set of all vectors Ci\v) 
with two more qubits up compared with the vector \v). 
If we by denote Kq the number of pairs of neighboring 
qubits both being down in the vector |i>) then the num- 
ber of elements in the set Mq is Kq. Equivalently, let 
us denote by M 2 to be the subset of all vectors Ci\v) 
with two more spins down as arc in the vector \v). The 
number of elements in the set M 2 is K 2 where K 2 is the 
number of pairs of neighboring qubits both being up in 
the vector \v) . Finally, let us denote by Mi the subset 
of all vectors C,|u) with equal number of qubits up as 
are in the vector \v). The number of elements set Mi is 
K\ and equals the number of pairs of neighboring qubits 
in the state \v) with exactly one qubit up. The following 
relations hold 

K 1 +K 2 + K 3 = N ; (D6) 
O.K0 + I.K1 + 2.K2 = 2K. 

We can rewrite in the following way Sj = s.ki where fcj 
is the relative sign of the vector Ci\v) according to the 
absolute sign of the vector \v) in the sum Eq. (|D1|) . 

Consider vectors belonging to the subset M . To find 
out the relative sign Sj between the vector Cj|u) and 
vector \v) we need to know the following distances 

1. The distance from any qubit up in the vector \v) to 
the j-th position: d(j, x) . 

2. The distance from the (j + l)-th qubit to any qubit 
up in the vector \v): d(j + l,x) . 

3. The distance of the j-th and the (j + l)-th qubit 
that is apparently one: d(j, j + 1) = 1 . 

The relative sign between the two vectors is then 

where summing over x means that we sum over all posi- 
tions of qubits up in the vector \v). 

Theorem 2 If the qubit in the vector \ v) that is equally 
distant from the j-th and the (j + l)-th qubit is down then 
the exponent in Eq. \DT\j is odd. 



Proof: There is only one qubit in the ring that is equally 
distant from the qubits on the j-th and the (j + l)-th po- 
sitions. If that qubit is down and taking any qubit in the 
vector \v) that is up and summing the distance from the 
j-th qubit to the given qubit and the distance from the 
(j + l)-th qubit to the same qubit we always get an odd 
number. As vector \v) contains even number of qubits 
up the sum is an even but to get the final expression in 
Eq. i|-L)7|) we have to add 1 therefore the exponent is odd. 

Theorem 3 If the qubit that is equally distant from 
the j-th and the (j + l)-th qubits is up in the vector \v) 
then the exponent in the sum in Eq. j-P7| ) is even. 

Proof: Follows from the previous statement. 

We know that we have Kq vectors in the set Mq. Denote 
ko to be the number of such states that they have the 
qubit which is equally distant from the corresponding j- 
th and (j + l)-th position up. Then we may say that the 
contribution of all vectors from the set Mq to the v-th 
component is 

- (ko - (K - k )).s , (D8) 

where the minus sign in front comes from Eq. (|D5ll . It is 
amazing that considering the other sets namely, Mi and 
M 2 , we have come to the same conclusion so that their 
contributions to the v-th component are 

- (kx - (Kt - kx)).s , (D9) 

from Mi and 

-(k 2 -(K 2 -k 2 )).s , (D10) 

from M 2 . Consequently, the v-th component of the vec- 
tor Hi\p) is a sum of Eq. Ip^l. Eq. |!23> and Eq. jPT0|) 
and reads 

(Kq + K X + K 2 - 2(k + h+ k 2 )).s . 

Now comes the crucial point. The following equation 
holds 

fc + ki + k 2 = K 

and together with Eq. (|D6|> the v-th component of Hj \p) 
is 

(H I \X)) v = (N-2K).s (Dll) 

c. We have proved that the left-hand side of 
Eq. HD3|) is equal to -(N - 2K).s and the right-hand 
side of the equation equals — (N — 2K).s. In other words 
the two expressions are equal for a given vector \v). Since 
we have not specified the vector \v), it holds for any vec- 
tor (see the discussion above) and we have proved that 
the X-state is an eigenstate of the Hamiltonian H with 
zero energy. 
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2. Density matrix of n neighboring qubits 

Theorem 4 Let N denote the total number of spins 
forming the chain and n be an integer. Further let N — 
2n + 1 and the system be in the X-state Eq. fDl\) . 

Then the density matrix of any sequence of n neigh- 
boring qubits is 



1 

PO - ^ In , 



where 1„ is the identity operator acting on the 2 n - 
dimensional Hilbert space of n qubits. 

The sequence of n neighboring qubits is a subset of all 
qubits in the chain such that only two "cuts" are needed 
to cut out the whole sequence from the chain. In what 
follows we will denote the set of n neighboring qubits by 
O. 

Proof of the Theorem 4 : We want to show that the 
density operator of n neighboring qubits is proportional 
to the identity operator acting on the 2™-dimensional 
Hilbert space. Denote by \{i, . . . k})o one of the basis 
state vectors of the system of n qubits from the set O 
where the set of indices {i, . . .} denote positions where 
the spins are up and all remaining spins are down. First 
we rewrite the X-state using this new basis as follows 

\X) = ]T {'.....,! o <>,,....,., o- (D13) 
{i,...k} 

where |df{i,...j})o is a state vector of the remaining n + 1 
qubits not belonging to the set O [41] and we sum over all 
sets of indices {«,... k} which means that we sum over 
all basis vectors of the system of n qubits. Then the 
Theorem 4 says that 

K 2 

o( a {k,...i}\a{i,...]})o = -^T s {k,...l},{i,...j} > (D14) 

where K is the norm of the X-state. 

In order to prove Eq. (|D14fl we need to know the form 
of the states \c*-{i,...j})o- From Eq. I|D1|) we have 

\X) = £ \{i,j,k,...}){-ljZ t >,W 

{i,j,-}e 

= E Kv--j}>° E i{ fc >-- J }>o 

{i,-j}e {k,...l}e 

+ E KV"J»o E K*>-l}>o 

{i,...j}o {k,...i\o 

where the sum over {i, . . .j} e means that we sum over 
all sets of indices with an even number of indices in each 
set, {j, . . . k}o means that we sum over all sets of indices 
with an odd number of indices in each set, du__n — 



J2i>j d (hj) where i,j S and d{i,...j},{fc,..j} = 

Sae{i j} be{k i}d{a,b). It follows that the states 
\ a {i,...j})o are 01 * ne following form. If the set {i, . . . j} 
contains an even number of indices then 

\<*{i,..j})o = E l{ fc '---0)o (D15) 

{k,...l}e 



X 



(_l) d { I ,..- 3 }+ d {fc..-.!}+ d { I ,..-j},{fc....l} 



(D12) while if the set contains an odd number of indices then 



l«{i,..,-}>6 = E K*."-l}>o ( D1 6) 

{k,...l}o 

X (— l) d {*."3}+ d {fc.--0+ d {i, -J}, {k,...l} 

Moreover, the norm of the X-state K can be easily cal- 
culated and the result is 



K 



™= E E i+ E E 1 

{i,...j}e [k,...l} e {i,...j}o{k,...l}o 



In what follows we consider three different possibilities. 

1. The two sets {i, . . .j}, {k, . . .1} are equal. If in the 
set {i, . . .j} is an even number of indices then 



oi a {i,...j} \ a {i,... j})o = E 
{*,...*}« 



1 = 2" 



else if in the set {i, . . . j} is an odd number of indices 
then 

o( a {i,-j}\ a {i,-j})d = E 1 = 2 " ■ 

{k,...l}o 

2. There is an even number of indices in the set 
{i, . . . j} while in the set {fc, . . . 1} the number of 
indices is odd. If the set {i, . . . j} contains an even 
number of indices from Ea. (|D15|l it follows that 
state \ocu...j\)o is a sum of basis state vectors with 
an even number of spins up. Further if the set 
{fc, . . . 1} contains an odd number of indices then 
from Eq. I|D16I) it follows that the state \a/k,...l})o 
is a sum of state vectors with an odd number of 
spins up and therefore the scalar product Eq. I|D14|) 
is zero i.e. 

o( a {i,-..j}\ a {k,...l})o = • 

3. The case that remains is when in the both sets 
{i, . . . j} and {fc, . . . 1} there is an even number of in- 
dices but the two sets are different 01 • The scalar 
product Eq. (|D14(1 reads as 

o( a {t...j}\a{k,...i})o = 

(_i)d{>,...j}+ d {™,-.-™}+ ti {i,...3}.{™,...™} 

{m, . . .n} e 

X (— l) rf {fc.---<}+ d {™----n}+ d {fc. ■••<}, {m,...n} 
= (— \) d {^--l}+ d {k,-l} 



x E (-D 

{m,...n} e 



2d, 



»}+<i{i,...i},{m,...n}+<i{fe,. ..!}.{">,...»} 
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Since we want to show that the scalar product of 
the two vectors is zero we drop the overall factor in 
front of the sum and using the relation (— l) 2k = 1 
we rewrite the last equation as follows 

o( a {i,...j}\ a {k,...l})o = 

^ (— ...n}+ d {k, ...l},{m,...»} 

{m,...n} e 

The problem of calculating the scalar product of 
two vectors has transformed into calculating the 
distances between the two sets. If an index a is in 
the both sets of indices {i, . . . j} and {k, . . . 1} then 
we can neglect it since we sum over distances from a 
to {to, . . . n} twice and the term (— l) 2d w,{™. ••■»} 
1 does not change the sign of the corresponding 
contributions. Therefore instead of calculating dis- 
tances '/;,...., ,„....„; and <2{fc,..j},{ m ,...„} , we cre- 
ate one set of indices 

D({i,...j},{k,...l}) = {,.....,! {/,'..../! 

- {i,...j}n{k,...l} , 

and then 

6i a {i,...j}\ a {k,-l})d = 

(_i) d D({i,... 3 },{fc,...i}),{ m , ...«} _ (di7) 

{m,...n} e 

It is important to note that the set 
D(\i, . . .j}, \k, . . always contains only even 
number of indices [42| • The last step is to calculate 
the distances. In order to calculate these distances 
we use the following strategy: We choose one of 
the positions from D({i, . . .j}, {k, . . . I}) such that 
it is the closest to the qubits not belonging to O 
and denote the position to be a. There are only 
two qubits that are equally distant from a and 
none of them belong to O. Denote their positions 
as b and c, where the position denoted as b is closer 
to the set O. Let us calculate the distances from 
D({i, . . . j}, {k, . . .1}) to b and c. If the distance is 

dD({i,...j},{k,...l}),{b} = C, 

then 

dD({i,...j},{k,...l}),{c} = C + (L - 1) , (D18) 

where L is the number of positions in 
D({i,...j},{k,...l}) as 

c) = d(i,b) + l; Vi e D({i, . . .j}{k, . . .1})/ a 
c) = d{i, c) ; i = a . 

Example: To make this clear let us consider a 
simple example N = 9, n — 4, O — {1,2,3,4} 
and D({i,...j},{k,...l}) = {2,4} see Fig. 00 The 



b=9 




a = 4 



FIG. 8: The ring of 9 qubits. O = {1, 2, 3, 4}. 

position that is the closest to the qubits not be- 
longing to O is 4 so that a = 4. The two posi- 
tions that are equally distant from 4 are 9 and 8. 
The position 9 is closer to O so that b = 9 and 
c = 8. Moreover the distances <i_D({v..j},{fc, ...;}),{&} 
and <2zj({i,...j},{fc,..j}),{c} are 

d D ({i,...j},{k,...i}),{b} = d(2, 9) + d(4, 9) = 1 + 4 , 
d D ({i,...j},{k,...l}),{c} = d(2, 8) + rf(4, 8) = 2 + 4 . 

The number of positions (indices) in 
D({i, . . .j}, {k, . . . I}) is equal to 2 and conse- 
quently the relation Eci. (|D18(l holds. 

Now we have everything necessary to calculate the 
scalar product Eq. I|D17|) . Choose one arbitrary set 
{m, . . . n} which contains neither b nor c. If there is 
an even number elements in the set {m, . . . n} then 
there are two contributions to the sum in Eq. I|D17|) 
namely 

<*£>({*,... j}, {fc,...l}),{m,...n} 

-f(-l) d D({ I ,...j},{fc,. ..!}), {b,c,m,...n} 

and these two have opposite signs [see Eq. QD18|l ]. 
It follows that these two contributions cancel each 
other. Equivalcntly, if there is an odd number of 
elements in the set {m, . . . n} then we have again 
two contributions 

(_l)<*D({<,...J},{*,...I}),{o.,m,...n} 
+ (— l) dj3 <{ i .- -..!}), {6,™.-™} 

with opposite signs and the two contributions are 
mutually cancelled. It follows that all contributions 
to the sum in Eq. (|D17|I mutually cancel and the 
result is zero. 
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To conclude, we have proved that 

o( a {k,...i}\®{i,...j})o = 2 ns {k,...i},{i,...j} , 

and it follows that the density operator of the n neigh- 
boring qubits is the identity operator up to a constant 
factor. 



3. Consequences 

d. Density matrix of a qubit Choose one of N qubits 
forming the Ising chain. Add n — 1 qubits to the chosen 
qubit such that the n qubits form a set of n neighboring 
qubits. Then according to Theorem 4 the density matrix 
of such system is 

1 

When we now trace over the n — 1 qubits that we added 
to the chosen qubit whose density matrix (state) we want 
to know, we obtain 

p = Tr„_ip„ = il , (D19) 



where 1 is the identity operator acting in the two- 
dimensional Hilbert space. 

e. Density matrix of a pair of qubits The derivation 
of the density operator of any two qubits follows the same 
steps as the derivation of the density operator of a single 
qubit. However, in this case we add only n — 2 qubits 
so that the set of n qubits consists only of n neighboring 
qubits 0|- According to Theorem 4 the density opera- 
tor of such system is p n = 1/2™ 1„. The corresponding 
density operator of the two qubits is obtained via tracing 
over the degrees of freedom of the n — 2 qubits that we 
added 



pij = Tr n - 2 pn = t I2 , (D20) 



where i, and j denotes the positions of the two a priori 
chosen qubits and I2 is the identity operator acting in 
the four-dimensional Hilbert space of the two qubits. 
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